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ABSTRACT 

Underwater  temperature  and  salinity  microstructure  can  lead  to  localized  changes  in  the  index  of  refraction 
and  can  be  a  limiting  factor  in  oceanic  environments.  This  optical  turbulence  can  affect  electro-optical  (EO) 
signal  transmissions  that  impact  various  applications,  from  diver  visibility  to  active  and  passive  remote  sensing. 
To  quantify  the  scope  of  the  impacts  from  turbulent  flows  on  EO  signal  transmission,  and  to  examine  and 
mitigate  turbulence  effects,  we  perform  experiments  in  a  controlled  turbulence  environment  allowing  the  variation 
of  turbulence  intensity.  This  controlled  turbulence  setup  is  implemented  at  the  Naval  Research  Laboratory 
Stennis  Space  Center  (NRLSSC).  Convective  turbulence  is  generated  in  a  classical  Rayleigh-Benard  tank  and 
the  turbulent  flow  is  quantified  using  a  state-of-the-art  suite  of  sensors  that  includes  high-resolution  Acoustic 
Doppler  Velocimeter  profilers  and  fast  thermistor  probes.  The  measurements  are  complemented  by  very  high- 
resolution  non- hydrostatic  numerical  simulations.  These  computational  fluid  dynamics  simulations  allow  for 
a  more  complete  characterization  of  the  convective  flow  in  the  laboratory  tank  than  would  be  provided  by 
measurements  alone.  Optical  image  degradation  in  the  tank  is  assessed  in  relation  to  turbulence  intensity.  The 
unique  approach  of  integrating  optical  techniques,  turbulence  measurements  and  numerical  simulations  helps 
advance  our  understanding  of  how  to  mitigate  the  effects  of  turbulence  impacts  on  underwater  optical  signal 
transmission,  as  well  as  of  the  use  of  optical  techniques  to  probe  oceanic  processes. 

Keyw^ords:  optical  turbulence,  laboratory  experiments,  ADV  profiler,  TKED,  TD 

1.  INTRODUCTION 

Turbulent  flows  are  an  intergral  part  of  the  natural  environment.  In  the  ocean,  the  mixing  that  accompanies 
turbulent  flows  is  an  important  part  of  the  global  energy  budget  and  of  localized  water  mass  exchanges.  Small- 
scale  temperature  and  salinity  variations,  i.e.  temperature  and  salinity  microstructure,  which  are  associated 
with  turbulent  mixing,  can  also  lead  to  localized  changes  in  the  index  of  refraction,  which  may  affect  underwater 
optical  properties.  This  phenomenon  has  been  demonstrated  to  have  the  potential  to  be  a  limiting  factor  in 
oceanic  environments,  affecting  optical  signal  transmissions  that  impact  various  naval  and  civilian  applications, 
from  diver  visibility  to  active  and  passive  remote  sensing.^ 

To  quantify  the  scope  of  the  impacts  from  turbulent  flows  on  optical  signal  transmission,  and  to  better  un¬ 
derstand  how  to  mitigate  these  impacts  in  the  underwater  environment,  we  implemented  a  controlled  turbulence 
setup,  where  convective  turbulence  is  generated  in  a  classical  Rayleigh-Benard  tank  (dimensions  are  5m  by  0.5m 
by  0.5m)  (Fig.  1).  The  goal  was  to  develop  a  setup  where  turbulence  strength  -  convective  turbulence  in  our 
case  -  can  be  adjusted,  while  also  being  characterized  in  terms  of  turbulent  kinetic  energy  dissipation  e  (TKED) 
and  temperature  variance  dissipation  rate  x  (TD).  This  controlled  turbulence  environment  can  then  serve  as  a 
testbed  for  optical  methods  to  mitigate  the  effect  of  turbulence  on  EO  signal  transmission. 
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Figure  1.  Laboratory  tank  at  NRLSSC.  The  photo  on  the  right  shows  the  Vectrino  profiler  ADV  and  the  CT  temperature 
probe. 


2.  METHODS 

Convective  turbulence  was  generated  in  the  RB  tank  by  heating  a  stainless  steel  plate  at  the  bottom  of  the  tank 
and  cooling  a  similar  plate  at  the  top  of  the  tank.  The  strength  of  the  convective  turbulence  in  the  tank  is  a 
function  of  the  temperature  difference  across  the  tank  and  can  be  characterized  in  terms  of  the  Rayleigh  number, 
defined  as 


^  gaAT 

Ra  =  - — - - 

i/Dt 


(1) 


Here,  g  is  the  acceleration  due  to  gravity,  a  is  the  thermal  expansion  coefficient,  AT  is  the  temperature 
difference  between  the  plates,  d  is  the  distance  between  the  plates,  u  is  the  kinematic  viscosity,  and  Dt  is 
the  thermal  diffusivity.  In  our  experiments,  Ra  ranges  from  1.5  •  10^®  to  around  4  ♦  10^^,  corresponding  to  a 
temperature  difference  across  the  plates  AT  of  6K  and  16i^,  respectively  (Table  1).  Note  that  this  temperature 
difference  is  the  difference  between  the  temperature  on  the  plates  proper,  the  temperature  gradient  that  sets  up 
in  the  interior  of  the  tank  is  much  less  pronounced. 

To  characterize  the  turbulent  environment  in  the  RB  tank,  the  turbulent  flow  is  quantified  by  high-resolution 
Acoustic  Doppler  Velocimeter  profilers  (Nortek  Vectrino  Profiler)  and  fast  thermistor  probes  (PME  high-resolution 
conductivity- temperature  (CT)  probe)  (Fig.  1,  right).  These  instruments  provide  high-resolution  velocity  and 
temperature  measurements,  at  100  and  6AHz  respectively.  The  ADV  and  CT  probes  were  mounted  in  the  RB 
tank  and  collected  time  series  of  high-resolution  velocity  and  temperature/conductivity  for  the  subsequent  esti¬ 
mation  of  TKED  and  TD  rates.  Data  were  collected  at  a  sampling  frequency  of  lOOHz  with  the  Vectrino  Profiler 
and  at  647^2;  with  the  CT  probe,  which  is  controlled  by  the  Nortek  Vector  ADV  in  our  setup.  A  high-speed 
imaging  camera  and  active  optical  target,  an  iPad  displaying  resolution  patterns,  were  placed  at  opposite  ends 
of  the  tank,  providing  an  optical  path  length  of  around  5m.  This  served  the  purpose  of  putting  the  turbulence 
data  into  the  context  of  measurements  of  optical  target  clarity. 

The  turbulent  kinetic  energy  dissipation  (TKED)  and  temperature  variance  dissipation  (TD)  rates  were  cal¬ 
culated  from  the  velocity  and  temperature  measurements  via  spectral  fitting  to  Kolmogorov  spectra  (for  velocity) 
and  Batchelor  spectra  (for  temperature)  and  compared  to  values  obtained  from  the  numerical  simulations  of  the 
RB  tank  for  comparable  Rayleigh  number  and  setup  (see  Section  2.1). 

2.0.1  Numerical  Tank 

The  laboratory  measurements  were  complemented  by  very  high-resolution,  non-hydrostatic  numerical  simula¬ 
tions.  These  simulations  provide  full  fields  of  temperature  and  velocity  for  the  estimation  of  turbulence  param¬ 
eters  and  their  impact  on  the  optics  (Fig.  2). 
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Figure  2.  Example  of  the  “Numerical  Tank”  used  to  simulate  Rayleigh- Benard  convection  and  emulating  the  laboratory 
tank  setup. 


The  numerical  experiments  (as  listed  in  Table  1,)  were  performed  with  the  open-source  Computational  Fluid 
Dynamics  (CFD)  package  OpenFOAM.  We  use  a  Large- Eddy  Simulation  (LES)  formulation,  with  the  traditional 
Smagorinsky  model  as  a  sub-grid  scale  model.  In  LES,  the  larger-scale  eddies  in  the  flow  are  explicitely  resolved, 
while  the  scales  smaller  than  the  grid-size  are  modeled.  We  run  at  two  resolutions:  a  centimeter  scale  simulation 
with  Ax  =  Ay  —  Az  =  1cm,  which  corresponds  to  1.25  million  grid  points  in  the  domain,  and  a  very  high- 
resolution,  millimeter-scale  simulation  with  Ax  =  Ay  —  5mm,  Az  =  2.5mm,  which  corresponds  to  20  million 
grid  points.  The  former  can  be  run  comfortably  on  a  modern,  high-performance,  dual  six-core  Linux  desktop, 
making  use  of  all  12  cores,  whereas  the  latter  requires  the  use  of  a  High-Performance  Computing  facility,  due  to 
the  high  computational  cost.  In  our  case,  the  code  was  run  on  96  processors  on  one  of  the  IBMs  at  the  Navy 
DoD  Supercomputing  Resource  Center. 

At  these  resolutions,  the  flow  is  very  well  resolved,  down  to  the  Kolmogorov  scales,  and  the  contribution  from 
the  sub-grid  scale  model  is  small,  especially  in  our  very  high-resolution  run.  In  this  paper,  we  show  results  from 
the  millimeter-scale  resolution  run,  unless  otherwise  noted. 

We  also  explored  cases  with  the  numerical  tank,  where  the  temperature  on  the  heating  (cooling)  plates  was 
not  uniform,  but  rather  decreased  (increased)  towards  the  center  of  the  tank.  The  plates  in  our  laboratory  setup 
are  heated  (cooled)  by  water  flowing  into  the  plates  and  these  experiments  helped  us  to  more  accurately  assess 
the  impact  of  heat  loss  over  the  plates  on  the  convective  turbulence  in  the  interior.  The  main  effect  appeared 
to  be  that  the  reduced  heat  input  lowered  the  Rayleigh  number  of  the  flow  and  thus  reduced  the  strength  of 
convection  in  the  tank.  To  help  quantify  this  result,  and  to  compare  simulations,  we  calculated  the  variance  of 
temperature  T  and  of  vertical  velocity  w  along  a  centerline  at  half-depth  in  the  tank.  These  values  are  shown  in 
Table  1. 

2.1  Calculation  of  TKED  and  TD 

2.1.1  Turbulent  Kinetic  Energy  Dissipation  rate  -  TKED 

Spectra  were  calculated  from  either  velocity  sections  (model)  or  from  velocity  time-series  (laboratory  data),  in 
which  case  Taylor’s  Frozen  Turbulence  hypothesis  (see,  f.  ex.,  Ref.  2),  k  =  had  to  be  applied  to  convert  the 
frequency  spectra  to  wavenumber  spectra.  Here,  k  is  the  wavenumber  (in  radm"^),  /  is  the  frequency  (in  s“^), 
and  U  is  the  mean  flow  past  the  sensor  (in  ms~^).  The  dominant  velocity  direction  is  vertical  in  our  case,  and 
we  use  the  mean  flow  associated  with  the  up/downdraft  past  the  sensor.  A  wavenumber  range  respresentative 
of  the  spectrum  in  the  ISR  is  chosen  and  the  TKED  rate  e  is  estimated  from 
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Exp#: 

1 

2 

3 

4 

5 

6 

7 

Ra 

4  •  10^^ 

4-10'^ 

2  •  10^^ 

1.5  -  10^^ 

4-101“ 

4  -  101“ 

4  -  lO^^i 

^Ttank 

16 

16 

8 

6 

16(14) 

16(12) 

16(8) 

^bott 

302 

302 

298 

297 

302 

302 

302 

Ttop 

286 

286 

290 

291 

286 

286 

286 

^'^plates 

no 

no 

no 

no 

±(1°) 

±(2°) 

±(4°) 

V8ir{w) 

3.6-10-5 

2.4-10-5 

1.1  -  10-5 

1.2- 10“^ 

2.1  - 10-5 

2.3-10-5 

1.9-10-5 

var(T) 

0.0060 

0.0040 

0.0009 

0.0007 

0.0024 

0.0014 

0.0016 

Resol. 

5a:2.5mm 

Icm 

Icm 

1cm 

Icm 

Icm 

Icm 

Table  1.  Table  of  numerical  experiments. 


Eii{ki)  =  ,  (2) 

where 

rOO 

Eii{ki)dki  =uf.  (3) 

Here,  Ea  is  the  one-dimensional  velocity  spectrum  of  component  i  [xi  =  a;  is  the  streamwise  component, 
^2  ?/  is  the  transverse  component,  and  =  2  is  the  vertical  component,  in  a  Cartesian  coordinate  system), 

and  uf  is  the  variance  of  the  signal.  C  is  a  constant  equal  to  ^  a  in  the  streamwise  direction  [x  and  fci)  and 
to  5  55  in  the  transverse  or  vertical  direction  {y,  2)  ;  a  «  1.5  is  the  Kolmogorov  constant  (see,  for  example, 3). 
Note  that  in  our  case,  since  we  use  the  vertical  velocity  as  the  mean  flow  U,  we  need  to  make  sure  to  use  the 
appropriate  C,  i.e.,  i|  a  for  the  vertical  direction  and  |  ^  q  for  the  horizontal  directions. 

2.1.2  Temperature  Dissipation  rate  -  TD 

For  the  estimation  of  TD  rates  from  the  temperature  sections  (model)  or  time  series  (laboratory),  we  calculated 
temperature  gradient  spectra  and  fit  them  to  the  theoretical  Batchelor  spectrum.^  This  technique  involves 
fitting  the  temperature  gradient  spectra  in  the  dissipative  high-frequency  range  of  the  spectrum,  which  is  poorly 
resolved  and  very  noisy  in  the  laboratory  data.  To  obtain  temperature  gradient  spectra  from  the  temperature 
time  series  in  the  laboratory,  we  again  invoke  Taylor’s  frozen  turbulence  hypothesis  and  use 

^ 

dx  ^  U~dt 

and  k  =  //{/  (note  that  in  this  case,  k  has  units  of  cpm)  .  For  isotropic  turbulence,  the  dissipation  rate  of 
the  temperature  variance  can  be  described  as 


dT  dT 

=  QDr  <fTAk)dk,  (5) 

where  Dt  is  the  thermal  diffusivity,  and  (k)  is  the  wavenumber  spectrum  of  the  alongpath  component 
of  the  temperature  gradient.^  To  ensure  that  U  varies  as  little  as  possible  over  the  segment  used  to  calculate 
temperature  gradient  spectra,  the  laboratory  spectral  estimates  were  done  from  Tsecond  data  segments  and  the 
results  for  x  were  averaged. 

We  explored  the  use  of  the  Maximum  Likelihood  Estimator  technique  (MLE)  of  spectral  fitting  to  our  data,^ 
which  works  very  well  for  the  wavenumber  spectra  directly  calculated  from  model  data,  without  the  use  of  the 
frozen  turbulence  hypothesis.  In  the  case  of  the  laboratory  data,  our  spectra  do  not  well  resolve  the  high- 
frequency  part  of  the  spectrum,  up  to  the  Batchelor  cut-off  wavenumber  and  we  instead  chose  to  integrate 
the  theoretical  spectrum,  according  to.^ 
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Figure  3.  Image  degradation  for  different  optical  turbulence  strength. 


2.2  Optical  Turbulence 

As  noted  previously,  the  turbulence  data  was  collected  in  the  context  of  optical  target  clarity  and  spatial  reso¬ 
lution.  High-definition  video  of  an  optical  chart  across  the  tank  was  taken  during  experiments  and  the  images 
assessed  for  image  degradation.  The  optical  chart  used  was  an  image  on  an  iPad  screen,  the  camera  was  a  Casio 
EXILIM  EX-Fl,  and  the  videos  were  shot  in  High-Definition  (HD),  with  30/ps  and  resolution  1920  by  1080 
pixels.  Note  that  the  ADVs  require  seeding  of  the  water  in  the  tank  to  collect  meaningful  velocity  data,  which 
affects  the  water  clarity.  Hence,  the  videos  shown  here  were  taken  immediately  preceeding  the  seeding  of  the 
water,  but  after  the  tank  had  reached  a  steady-state  with  respect  to  convective  turbulence. 

Images  from  the  laboratory  show  the  image  degradation  due  to  optical  turbulence  for  a  case  of  strong  (Fig. 
3,  center,  AT  ^  6K)  and  extreme  optical  turbulence  (Fig.  3,  right,  AT  ~  16K).  Here,  particle  scattering  is 
secondary  to  the  changes  in  the  index  of  refraction  (lOR)  due  to  temperature  microstructure.  Note  that  the 
effect  of  optical  turbulence  is  more  pronounced  at  the  higher  spatial  frequencies.  To  quantify  the  extent  of  image 
degradation  from  optical  turbulence,  we  apply  the  Structural  Similarity  Index  Method  (SSIM).®  The  SSIM  is 
an  image  quality  metric  and  measures  the  similarity  between  two  images,  when  one  is  considered  the  “perfect 
quality”  reference  image  (Fig.  3,  left). 


3.  RESULTS 

3.1  Image  degradation 

Results  from  applying  the  SSIM  metric  to  the  images  collected  at  the  different  turbulence  strengths  shown  in 
Figure  3,  show  increased  degradation  for  stronger  optical  turbulence,  i.e.,  higher  Rayleigh  number  (Fig.  4)  . 

3.2  Model  Temperature  Variations  and  Circulation 

To  characterize  the  turbulence  in  the  tank  that  is  associated  with  the  image  degradation  shown  in  Figure  4, 
we  analyze  the  results  from  the  laboratory  measurements  and  the  numerical  tank.  In  the  case  of  the  numerical 
experiments,  we  want  to  make  sure  that  the  model  has  reached  a  steady-state  and  that  the  convective  flow  is 
fully  developed.  To  investigate  this  time  evolution  of  the  convective  turbulence  in  the  domain,  we  look  at  a 
time  series  of  the  average  value  of  the  variance  along  a  spatial  section  taken  from  the  center  line  of  the  tank, 
for  both  temperature  and  vertical  velocity.  The  results  show  that  the  model  reaches  a  steady-state  after  about 
15  —  20  minutes  (Fig.  5,  data  shown  are  from  the  cm-scale  resolution  run).  The  variability  after  this  time  can 
be  attributed  to  the  dynamics  of  convective  cells  and  the  circulation  in  the  tank. 

These  convective  cells  are  also  apparent  in  temperature  and  velocity  fields  from  the  model  (along  the  section 
shown  in  Fig.  2)  and  show  the  convective  cells  that  develop  in  the  tank  (Fig.  6).  The  size  of  these  convective 
cells  is  governed  by  the  height  of  the  tank,  and  we  count  about  four  to  five  up-  and  downdrafts.  The  index  of 
refraction  (lOR)  in  the  tank  can  be  calculated  from  the  temperature  fields'^  and  also  shows  variability  associated 
with  the  convective  cells  in  the  tank.  Successive  cross-sections  of  temperature  taken  halfway  along  the  tank  (at 
X  =  2.5m)  also  illustrate  the  circulation  (Fig.  7)  and  variability  associated  with  this  dynamic  system.  While 
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Figure  4.  Image  degradation  for  different  (optical)  turbulence  strength  as  shown  in  Fig.  3  quantified  by  applying  the 
SSIM  method. 


Proc.  of  SPIE  Vol.  9111  911113-6 


Downloaded  From:  http://proceedings.spiedigitalIibrary.org/  on  06/04/2014  Terms  of  Use:  http://spiedl.org/ternis 


0.5  pr- 

£ 

w 

't  ^ 

•  -A. 

m 

£  0.25 
o 

....»-*  . . . 

..  . 

^*293 

0  1  2  3  4  5 

along-tank  distance  (m) 


_E 

£ 


2  3 

along-tank  distance  (m) 


't  V 

-T - — 

■”  T  "  -1  V - 

— ^ - y-' . V . . 

A  :r^»  ' 

1  1.3352 
1.335 

I  i  ...  :  , 

1.3348 

M.3346 

1 

0 

1 

2 

3 

4 

< 

along-tank  distance  (m) 


Figure  6.  Model  results  at  time  t  =  1050s,  showing  a  sideview  of  temperature  T  (top),  vertical  velocity  (Uz)  (middle)  and 
the  Index  of  Refraction  (lOR)  calculated  from  the  temperature  field  (bottom).  The  data  is  shown  on  the  center  plane  of 
the  tank. 


Figure  7.  Snap  shots  of  a  tank  cross  section  of  temperature  (in  °C)  illustrate  the  variability  associated  with  the  convective 
cells  and  the  dynamic  environment  in  the  tank  (x  =  2.5m;  at  times  t  —  850s,  875s  and  950s,  from  left  to  right). 


the  velocities  associated  with  the  flow  are  small  (0(2cm/s),  the  overturning  still  appears  vigorous  when  we  look 
at  the  temperature  microstructure.  The  cross-sectional  circulation  is  also  apparent  when  looking  at  streamlines 
from  the  simulation.  Figure  8  shows  the  updraft  on  one  side  of  the  tank  near  the  center  and  downdrafts  on  the 
opposing  side. 

We  can  again  calculate  the  lOR  from  the  model  temperature  fields,  and  the  cross-sections  shown  in  Fig.  9 
(top)  mirror  the  variability  seen  in  the  temperature.  To  assess  the  effect  on  an  optical  beam  across  the  tank, 
which  experiences  the  cumulative  effect  of  the  turbulent  variations  in  its  path,  we  averaged  the  lOR  at  each 
point  in  the  y-z  plane  over  the  entire  length  of  the  tank.  As  can  be  expected  from  the  averaging  procedure, 
this  notably  reduces  the  variability  seen  in  the  cross-sections.  Since  the  beam  passes  through  different  layers 
of  turbulence  along  its  path,  each  of  which  varies  in  time  and  space,  averaging  of  the  fields  may  not  be  the 
appropriate  approach  to  estimate  the  extent  of  beam  distortion  due  to  cumulative  changes  in  lOR.  Rather,  a 
more  involved  approach  directly  tracing  the  beam’s  path  through  the  medium  may  need  to  be  adopted  in  future 
work. 


Proc.  ofSPIE  Vol.  9111  911113-7 


Downloaded  Prom:  http://proceedings.spiedigitallibrar\.org/  on  06/04/2014  Terms  of  Use:  http://spiedl.org/terms 


uz 


-0.02455  0.02274 


UZ 


-0.02455  ’  0.02274 


Figure  8.  Model  results  at  time  t  =  1050s,  showing  streamlines  in  the  domain  colored  by  vertical  velocity  Uz  (top).  The 
bottom  plot  shows  the  same  but  looking  at  it  from  the  other  side  of  the  tank  (x-  and  y-axis  reversed).  The  updrafts  and 
downdrafts  illustrate  that  there  is  a  circulation  in  the  tank  in  the  cross-sectional  direction. 
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Figure  9.  Top:  The  lOR  calculated  from  the  snap  shots  of  the  tank  cross  section  of  temperature  shown  in  Fig.  7  illustrate 
the  variability  in  terms  of  lOR.  Bottom:  Since  an  optical  beam  across  the  tank  experiences  the  cumulative  effect  of  the 
turbulent  variations  in  its  path,  we  averaged  the  lOR  at  each  point  in  the  y-z  plane  over  the  entire  length  of  the  tank. 
As  can  be  expected,  the  averaging  notably  reduces  the  variability  seen  (at  times  t  =  850s,  875s  and  950s,  from  left  to 
right). 
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Figure  10.  Energy  spectra. 
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Figure  11.  Temperature  gradient  spectra. 

3.3  Spectra  -  Laboratory  and  Model 

As  described  in  Section  2.1,  we  calculate  TKED  and  TD  rates  from  both  laboratory  experiment  and  model 
to  characterize  the  convective  turbulence  in  the  tank.  Energy  spectra  E{k)  comparing  lab  (Fig.  10,  left)  and 
model  (Fig.  10,  right)  show  that  the  numerical  model  (especially  at  mm-scale  resolution)  adequately  captures 
the  turbulence  spectrum  in  the  inertial  subrange  and  allows  for  the  estimation  of  the  TKED  rates  e,  which  are 
of  0{lO~'^W/kg)  for  the  experiments  shown. 

The  energy  spectra  do  not  show  a  significant  difference  in  TKED  between  the  experiments  with  different 
turbulence  strength.  As  expected,  the  model  energy  spectra  show  a  drop-off  at  the  higher  wavenumbers  towards 
the  dissipative  range,  which  is  resolution  dependent. 

While  the  model  energy  spectra  accurately  capture  the  inertial  subrange,  at  levels  comparable  to  the  labo¬ 
ratory,  the  temperature  gradient  spectra  show  that  the  model  (Fig,  11,  right,  as  well  as  faint  blue  on  left)  has 
more  difficulty  reproducing  the  temperature  gradients  correctly.  It  is  these  temperature  gradients,  which  are 
predominantly  needed  to  assess  optical  turbulence  The  spectra  illustrate  that  for  different  optical  turbulence 
strengths,  while  the  mechanical  turbulence  is  comparable,  the  temperature  gradient  fields  show  a  marked  differ¬ 
ence.  Optical  turbulence  is  mainly  due  to  temperature  variations  affecting  the  lOR  of  water,  and  to  adequately 
describe  the  effect  on  the  optics,  it  is  critical  to  resolve  the  temperature  gradients.  The  value  of  x  (which  was 
calculated  following  Ref.  4)  is  also  strongly  dependent  on  model  resolution.  The  contribution  of  the  sub-grid 
scale  temperature  fluctuations  needs  to  be  estimated  to  directly  compare  the  temperature  gradient  spectra  be¬ 
tween  model  and  laboratory.  Estimating  the  model  sub-grid  scale  contribution  to  x  eliminating  potential 
errors  associated  with  the  use  of  Taylor’s  Frozen  Turbulence  hypothesis  for  the  calculation  of  x  (and  e)  from  the 
laboratory  data,  are  the  subject  of  ongoing  work. 

4.  SUMMARY  AND  DISCUSSION 

In  order  to  study  the  impact  of  temperature  microstructure  on  underwater  optical  signal  transmission,  we  per¬ 
formed  experiments  in  a  controlled  laboratory  environment  complemented  by  high-resolution,  non-hydrostatic 
numerical  simulations.  The  goal  was  to  develop  a  setup  where  turbulence  levels  can  be  controlled  and  fully  char¬ 
acterized.  This  setup,  which  allows  for  repeatable  experiments  under  controlled  conditions,  can  help  understand 
processes  involved  in  optical  turbulence  and  provide  a  platform  for  the  testing  of  optical  techniques  to  mitigate 
turbulence  effects  underwater.  Optical  turbulence  is  mainly  due  to  temperature  (or  salinity)  variations  affecting 
the  lOR  of  water,  and  to  adequately  describe  the  effect  on  the  optics,  it  is  particularly  critical  to  resolve  the 
temperature  gradients.  This  can  present  a  challenge  in  both  the  laboratory  and  the  model,  due  to  noise  and  res¬ 
olution  requirements,  respectively.  Further  work  is  needed  to  adress  questions  related  to  sub-grid  contributions 
in  LES  to  the  rate  of  temperature  variance  dissipation. 
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